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ABSTRACT 


This  report  is  a  continuation  of  the  previous  study  on  simulation 
parameters.  Additional  cases  studied  in  this  report  indicate  that  the  minimum 
number  of  points  per  wavelength  or  per  correlation  length  whichever  is  the 
smaller  is  around  5,  the  minimum  number  of  surface  samples  is  25  and  the 
minimum  width  of  the  illuminated  area  is  7  X.  or  7L  whichever  is  the  larger. 
These  studies  are  still  incomplete  in  that  we  have  not  considered  enough 
variations  in  the  surface  rms  height,  surface  rms  slope  and  the  ratio  of  the 
incident  wavelength  to  the  surface  correlation. 
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1.0  INTRODUCTION 


This  is  the  second  quarterly  report  on  the  Contract  N00014-87-K-0751 . 
In  our  first  report  [Fung  et  al,  1987]  we  have  examined  the  problem  of  surface 
generation,  the  calculation  of  the  backscattering  coefficient  and  the  conditions 
on  the  selection  of  simulation  parameters  for  certain  types  of  surface.  This 
quarter  we  have  an  update  on  our  surface  simulation  program  and  we  continue 
our  studies  on  the  selection  of  simulation  parameters  for  additional  surface 
types.  This  means  that  finer  divisions  of  cases  are  considered  to  realize  greater 
efficiency  in  the  simulation  program.  The  need  for  an  efficient  program  was 
pointed  out  in  our  previous  quarterly  report. 

In  Appendix  A  we  provide  an  updated  version  of  our  surface  generation 
program  where  a  more  efficient  fast  Fourier  transform  algorithm  is  used  for 
filter  weight  calculations  and  surface  generation. 

Based  upon  our  studies  in  the  previous  quarter  we  have  found  that  we 
need  a  minimum  of  4-5  points  per  wavelength  X  ;  with  a  Gaussian  antenna 
pattern  exp(-x:/g2)  we  need  an  effective  illuminated  area  2g  equals  to  about  3 
correlation  lengths  L  and  the  edge  of  the  illuminated  area  should  be  taken  to  be 
the  point  where  the  antenna  pattern  has  decayed  to  10 3.  In  all  cases  50  surface 
samples  were  used  to  obtain  a  statistically  averaged  value.  In  these  studies 
surface  parameters  with  rms  surface  height  o  of  0.5  unit  and  X  =  L  =  2  units 
were  used.  In  this  quarter  we  continue  the  same  study  by  examining  the  effects 
of  surface  parameters  on  the  selection  of  D,  the  width  of  the  illuminated  area, 
and  the  number  of  points  per  wavelength.  X ,  or  correlation  length,  L, 
whichever  is  the  smaller  and  the  total  number  of  samples.  N,  under  different 


combinations  of  the  following  surface  conditions:  (1)  L>  X  ,  (2)  L<  X,  (3) 
normalized  rms  height,  ks  =  2tt  c/  X  <1  ,(4)  ka  >  1.  The  same  antenna 
pattern  is  used  in  these  studies  as  in  the  previous  quarter. 
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2.0  SELECTION  OF  SIMULATION  PARAMETERS 


In  what  follows  N  the  number  of  surface  samples,  D  the  width  of  the 
illuminated  area,  and  the  number  of  points  per  correL.on  length  L  or  per 
wavelength  X  whichever  is  the  smaller  are  considered  under  different  surface 
conditions.  It  is  anticipated  that  D  is  a  function  of  X  and  L  but  it  is  not 
immediately  obvious  whether  it  should  depend  only  on  X  or  L  whichever  is 
the  larger.  N  is  expected  to  be  a  function  of  the  normalized  rms  height  kc. 
The  number  of  points  per  wavelength  or  correlation  length  whichever  is  the 
smaller  may  depend  on  the  relative  size  of  X  and  L.  It  is  possible  that  the 
required  number  of  points  per  correlation  length  may  be  different  from  the 
number  of  points  per  wavelength  . 

2.1  Number  of  Surface  Points  Per  Wavelength 

To  investigate  the  number  of  points  per  wavelength  we  must  make  sure 
that  the  requirements  on  other  surface  parameters  are  satisfied.  In  Figure  1 
we  have  chosen  L  greater  than  X  ,  D  greater  than  10  L  and  an  N  equal  to  50. 
We  shall  show  in  subsequent  sections  that  both  D  and  N  are  sufficiently  large. 
In  addition,  since  L  is  greater  than  X,  it  automatically  contains  enough  number 
of  points  if  X  contains  enough  number  of  points.  The  number  of  points  per 
wavelength  considered  ranges  from  1.3  to  8  in  Figure  1.  From  Figure  1  (a)  it 
is  seen  that  the  backscattering  coefficient  curves  begin  to  converge  when  the 
number  per  wavelength  exceeds  2.6.  The  difference  between  4  points  and  8 
points  per  wavelength  is  around  2  dB  or  less.  The  same  observation  can  be 
made  when  the  rms  surface  height  is  increased  from  0.3  to  0.6  in  Figure  1(b). 
Further  increase  in  the  rms  surface  height  o  to  1.2  in  Figure  1(c)  does  not 
seem  to  cause  any  change  in  the  number  of  points  required  per  wavelength. 
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Thus,  we  conclude  that  depending  upon  the  amount  of  error  one  can  tolerate 
the  required  number  of  points  per  wavelength  should  be  between  4  and  8. 

2.2  Number  of  Surface  Points  Per  Correlation  Length 

To  investigate  the  number  of  surface  points  per  correlation  length  we 
choose  the  wavelength  X  greater  than  the  correlation  length  L.  In  Figure  2(a) 
the  number  of  surface  points  per  correlation  length  is  varied  from  0.8  to  8.  It 
is  seen  that  convergence  in  the  backscattering  curves  occurs  when  there  are 
more  than  two  points  per  correlation  length.  It  appears  that  4  to  8  points  per 
correlation  length  seem  to  produce  satisfactory  results.  Then  in  Figure  2(b) 
the  rms  surface  height  is  increased  from  0.15  to  0.3  unit.  Similar  observations 
remain  applicable.  Next,  we  keep  other  parameters  the  same  and  increase  the 
correlation  length  from  1  to  2  units  and  maintain  the  same  surface  slope  by 
doubling  the  surface  rms  height  in  Figure  2(c)  and  2(d).  In  these  cases 
comparable  accuracy  to  previous  cases  can  be  obtained  by  requiring  more  than 
5  points  per  correlation  length.  It  is  not  clear  whether  this  somewhat  more 
stringent  requirement  comes  from  a  larger  L  or  a  larger  rms  height  or  the 
higher  order  surface  frequency  components  resulting  from  the  use  of  finite 
record  length  in  surface  generation. 

2.3  Width  of  the  Illuminated  Area 

Two  different  cases  are  considered  when  study  the  width  of  the 
illuminated  area  :  (1)  when  L  >  X  and  (2)  when  X  >  L.  In  case  (1)  we  want  to 
examine  how  many  L  should  be  contained  in  the  illuminated  area  D  so  that  the 
scattering  coefficient  calculated  agrees  with  that  from  an  infinitely  large 
surface.  Similarly,  in  case  (2)  we  want  to  find  out  how  many  X  should  be 
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contained  in  D  to  achieve  the  scattering  coefficient  from  an  infinite  area. 

Figure  3  shows  the  case  with  L  >  X  .  All  surface  parameters  are  given  in 
the  figure  where  K  is  the  wave  number  and  'sig'  stands  for  a  .  The  width  of  the 
illuminated  area  is  varied  from  D  =  5L  to  8L.  It  is  seen  that  the  backscattering 
curves  appear  to  converge  when  D  is  greater  than  7L.  In  all  the  computations 
8  or  more  points  per  wavelength  are  used  and  a  minimum  of  45  scattering 
coefficient  samples  has  been  averaged.  The  minimum  required  size  of  the 
illuminated  area  seems  to  be  7L  depending  upon  the  amount  of  error  one  can 
tolerate.  In  addition,  the  theoretical  scattering  coefficient  predicted  by  the 
Kirchhoff  surface  scattering  model  [Ulaby  et  al,  1982]  is  also  plotted  in  Figure 
3  to  provide  a  reference.  The  polarization  used  in  this  report  is  HH 
polarization.  Since  the  Kirchhoff  scattering  model  has  no  polarization 
dependence,  in  general  we  expect  fair  agreement  for  small  angles  of  incidence 
up  to  about  25  degrees  and  higher  theoretical  than  simulated  values  at  large 
angles  of  incidence  [Fung  and  Pan,  1987], 

Figure  4  shows  the  case  with  X  >L.  All  surface  parameters  are  given  in 
the  figure.  Since  X  >  L  we  seek  the  width  of  the  illuminated  area  D  in  terms  of 
X  and  vary  D  from  D=5  X  to  9  X.  When  D  is  greater  than  or  equal  to  7  X  , 
the  computed  scattering  coefficients  are  practically  the  same  indicating  that 
convergence  of  results  has  occurred.  Thus,  the  minimum  width  of  the 
illuminated  area  required  in  this  case  is  7  X  .  Note  again  that  we  have  used  N  = 
45  which  is  larger  than  the  minimum  required  values  for  the  surface  samples 
(to  be  demonstrated  in  the  nest  section)  and  a  minimum  of  8  points  per 
correlation  length  in  all  the  calculations. 


2.4  Number  of  Surface  Samples  Averaged 


In  Figure  5  we  illustrate  the  requirement  on  the  number  of  surface 
samples  needed  to  achieve  an  acceptable  averaged  value  for  the  backscattering 
coefficient.  Because  the  selected  case  is  a  slightly  rough  surface  with  a  small 
normalized  rms  surface  height,  ko  of  0.24,  there  is  a  significant  contribution 
from  the  coherent  component  of  the  backscattering  coefficient.  This  explains 
the  peaking  near  the  nadir  region.  For  the  case  shown  convergence  appears  to 
occur  after  the  number  of  samples  N  reaches  25.  Here  again  to  be  sure  that 
other  surface  and  system  parameters  do  not  affect  the  results  obtained  we  use 
values  which  exceed  the  minimum  requirements  for  all  of  them.  These  values 
are  given  in  the  figure  except  the  number  of  points  per  correlation  length 
which  is  8. 

In  Figure  5  we  have  plotted  the  total  scattering  coefficient  which  is  the 
sum  of  the  coherent  plus  the  incoherent  scattering  coefficients.  When  kct  is 
larger  than  unit)  ihc  coherent  contribution  becomes  negligible.  For  this  case 
and  the  case  of  L>  X  we  illustrate  the  requirement  on  N  in  Figure  6.  Here  N  is 
varied  from  5  to  30.  It  is  seen  that  after  N  exceeds  20  a  gradual  convergence  is 
obvious.  Depending  upon  the  amourii  of  error  one  can  accent  N  -  25  seems  to 
be  the  minimum  number  of  samples  needed. 
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3.0  DISCUSSIONS 


From  the  above  studies  we  note  that  the  minimum  number  of  points  per 
wavelength  or  per  correlation  length  whichever  is  the  smaller  is  around  5.  the 
minimum  number  of  surface  samples  is  25  and  the  minimum  width  of  the 
illuminated  area  is  7  X  or  7L  whichever  is  the  larger.  The  minimum  number 
of  points  per  X  or  L  may  have  a  dependence  on  the  ratio  of  X/L.  This  has  not 
been  studied  extensively.  In  the  study  of  the  width  of  the  illuminated  area  we 
also  have  not  considered  enough  variations  in  surface  roughness  that  is  a 
change  in  ka  and  g/L.  In  addition,  there  may  be  a  dependence  on  how  perfect 
is  the  generated  surface  meeting  its  statistical  description.  For  example,  the 
difference  between  the  required  number  of  points  per  correlation  length  in 
Figure  2(b)  and  2(d)  may  be  caused  by  possible  presence  of  high  frequency 
components  in  the  surface  not  reflected  by  the  correlation  length  parameter. 
To  better  understand  this  problem  some  smoothing  operation  should  be 
performed  on  the  surface  to  ensure  that  there  are  no  unwanted  high  frequency 
surface  components. 

For  a  statistically  described  surface  there  is  always  a  finite  probability 
for  it  to  take  on  values  substantially  different  from  its  mean  even  when  the 
surface  statistics  perfectly  match  its  required  specification..  Furtheimore. 
since  we  can  only  deal  with  finite  samples,  the  averaged  result  may  be 
affected  substantially  by  one  or  two  large  samples  with  very  low  but  finite 
probability  of  occurrence.  Thus,  special  care  must  be  exercised  in  estimating 
the  mean  scattering  coefficient  especially  at  large  angles  of  incidence. 
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5.0  FIGURE  LEGENDS 


1.  Dependence  of  the  hackscattering  coefficient  on  the  number  of  points  per 
wavelength  when  N  is  50  .  the  rms  surface  height  is  (a)  0.3  unit  (b)  0.6  unit 
and  (c)  1 .2  units. 

2.  Dependence  of  the  backscattering  coefficient  on  the  number  of  points  per 
correlation  length  L  with  N  =  50  and  when  (a)  the  rms  surface  height  is  0.15 
unit  and  the  correlation  length  is  1  unit  (b)  the  rms  surface  height  is  0.3  unit 
and  the  correlation  length  is  1  unit  (c)  the  rms  surface  height  is  0.3  unit  and  the 
correlation  length  is  2  units  and  (d)  the  rms  surface  height  is  0.6  unit  and  the 
correlation  length  is  2  units. 

3.  Dependence  of  the  backscattering  coefficient  on  the  width  of  the  illuminated 
area  D  when  the  correlation  length  is  larger  than  the  electromagnetic 
wavelength. 

4.  Dependence  of  the  backscattering  coefficient  on  the  width  of  the  illuminated 
area  D  when  the  correlation  length  is  smaller  than  the  electromagnetic 
wavelength. 

5.  Dependence  of  the  backscattering  coefficient  on  the  number  of  surface 
samples  N  when  the  correlation  length  is  smaller  than  the  electromagnetic 
wavelength. 

6.  Dependence  of  the  backscattering  coefficient  on  the  number  of  surface 
samples  N  when  the  correlation  length  is  larger  than  the  electromagnetic 
wavelength. 
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igUre  1  (a)  Dependence  of  the  oacs scatter ing  coefficient  on  the 
nuracer  of  points  per  wavelength  when  N  is  30,  tne 
rms  surface  height  is  0.3  unit 
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Dependence  of  the  oack scatter ing  coefficient  on  tne 
numoer  of  points  per  wavelength  when  N  is  L>0,  the 
runs  surface  height  is  1.2  units 
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Figure  2(a)  Dependence  of  the  oacksca tter ing  coefficient  on  tne 

nurmaer  of  points  ^er  correlation  lengtn  L  with  L  =  -<0 
and  when  the  rras  surface  height  is  0.15  unit  and  the 
correlation  length  is  1  unit 
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Figure  2(d)  Dependence  of  the  back  scatter ing  coefficient  on  the 

number  of  points  per  correlation  length  L  with  N  =  ^0 
and  when  the  rums  surface  height  is  0.3  unit  and  the 
correlation  length  is  1  unit 
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Figure  2(c)  Dependence  of  the  oack sea tter ing  coif f icient  on  tne 
number  ol  points  per  correlation  length  L  with  N  = 
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correlation  length  is  2  units 
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Figure  4  Dependence  of  the  Dack scatter ing  coefficient  on 
the  width  of  the  illuminated  area  D  when  the 
correlation  length  is  smaller  than  the  electro¬ 
magnetic  wavelength 
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Figure  6  Dependence  of  the  backscattering  coefficient  on  the 
numoer  of  surface  samples  N  when  the  correlation 
length  is  larger  than  the  electromagnetic  wavelength 
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V .  at-.. 

V  v  •»  *  w 

•>L.  ^  '»•  \ 


-  V  '  , 


\  "■  A  "-£  rOUP»F&  TRANSFORM  THt  .•  N^^NCt>ON  ; 
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IP 


M  '  =  SU  J) 

~  SORT  (,  S'.  I  .v1 ' ' 

K  *"  ^  Z 


NOS'  ' 


ONT 


N  r 


NT  ‘.  -~ 


'  T 


_a:-  ;•  :'an.;e  r:  ."■• 

Hr  DC*^:NANT  co:*  >, 

-0: l.  r  POI"  -OK-  ■  ;  '■' ! 


A>  \  A  '  Si  A  • 


=  1  .NWYh 

A-  \  A  >  -  ,NWYH-  :  *  '  )  =  $(  !  ,1 ) 

A  ■  \  A  >  “  ,NWYH  *  I  -  1  )  =  S(  I  ,1) 

o; 

j,';  1  4  '  .=  ’  NWVH 
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► 
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WkiNWAH-l-  !  .NWYH-J*  !  )=S(I,J) 

'  5C  CONTINUE 
160  CONTINUE 
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NNX=NN*N> 
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>  t  M  2  —  ^ ;  N 1  ' 
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END 
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CALL  FCJ’-.  -NT  ,5>  1  '•  f  •  -NT*  :  . . 
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r- 
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C 
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call  r:P~:  m.e' 
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?;2  =  2  *P'I 
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: - ;  =  ■  n 
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C 
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NN=4*N 

20  50  ■*  =  1  ,N4P0w 
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v 
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50 
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2  ■  ON  50( 2),8  K2) 

DO  1  0  ►  =  T  INT 


RETURN 

END 
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\  V-  • 

R‘\i  LO  ,L(  I  C 1 ),( L 
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DO  40  K  =2, 15 

!Fu(K-n-:'  io,:o,30 

10  Li  * '  i  i=2 

20  un;=2 

00  T0  40 

50  L,kl  =  L\N-i^2 

40  CONT|NUE 
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DC  '  2  C  i ,„  ~  l  : 

DO  !  20  j?  =  .  7  L?  ,l  7 
DO  '20  ,:4=j6, l3,l6 
DO  ’  20  v  i  0=09, L 1 0,19 
DO  ’20.;  1  =J  1 0,L  I  1  ,L  10 
DC  120  j ' 2=J1  1  ,L12,L!  1 
DO  ’ 20 O’ 3=01 2, L  13.L  12 
DO  1 20  o '  4=..., ;  3  i  1 4  ,l  i  3 
DO  1 20  ;thE”  =  '  -i  w  1  c  '4 
Rw2=.RHET-2 


DC  60  f  -  .;n; 

T0=B0^ -*P2(>  ' 

71  =BUk  '■♦63v  *■ 
B2(K)=BOOO-B2(K) 
B5viO=Bl(K)-B3(lO 
?(h>  )=T0  +  T  1 
8  l  (  k)=T0-T 1 
CONT'NUF 


F ,  NN-4 )  120,120.  0 
K0=iNT*4+  1 
_=<>;♦;  fsiT  -  ■ 

■C  80  k  =  6C  u 

PP=P7*‘  p  '  ( »  i - g; 3 1  k"  >  > 

p  ,  a  ►  ~  *  ;  n  1  .  k  i  ♦  C’  5  \  k  )  ) 

P3>>  ;=?:  k!  +  f 1 
g  i ;  k  )  =  P ;  -  p  2'  *  1 

52( f  )=50> ‘  '--Of 

C  V  ,  '  =  P  ;  ¥  '  *  P  C- 

..On^'NP- 


APO=THZ*PlCVN 

:  ■  =coSv  a-c-n 
'  =  8 '  s> aP.-  ■ 

•  ->  =  '  1  t  **2 


'?=C1*8 ’ *  62*C ' 

Ni*  4  =  "■,:*  4 
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0'=5i(o;*S!+B5^.)*Cl 
T2=82(J)*C2-B6(K  )*52 
T,t'  =  p2(J)*52  +  B6(tr  )*C2 
’  3 = c  3(  J  )*C3  -  5  7(  K  ;*S3 

’-.C;  f  ;'4C7*0*?U']*[7 

=  p  i":  i  ♦  7  j 

~  4  =  R  -*•  k  .  ♦  T  f. 


T i =r > ♦T3 
TS=a5+T7 
T  3 = k ' -13 
T?=P5-T7 
BOt  v'  )=TO+  T  I 
B7i k.  )=T4+T5 

5  i  )=j q-t  i 
p  :u  i=TC-T4 

f  ■  3 1  K  I  =  T  KV  ♦  T  y, 

5  -R  K  '■ = T  Z  * " ~ 
P3u')=T3-Tt 
'•CO  CONTINUE 
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v  =j'-: 
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110  JI=2*jP - i 
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!  20  CONTINUE 
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L^l'v-'  'LS.L  -0  ■).!  L5,L(  ’  ’  .".iL4.Lv  ’2' 
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4C  CONT-NME 


p:OV\=R 
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53=C2*S1  ♦  :2*C ■ 
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*  ?=v|!*INT4+  1 
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55v k ;=v ♦! 
£2\  A  =  (  T0--4'*C2-C 
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e:.k 
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